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Abstract 

This is the first of two papers in which we study real and virtual photon-proton 
scattering in a nonperturbative framework. We classify different contributions 
to this process and identify the leading contributions at high energies. We then 
study the renormalisation of the photon-quark-antiquark vertex that occurs in the 
leading contributions. We find something like the dipole picture in one of these 
contributions but also find two correction terms which can potentially become large 
at small photon virtualities. In the second paper we will discuss the additional 
approximations and assumptions that are necessary to obtain the dipole model of 
high energy scattering from the results found here. 
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1 Introduction 



The dipole picture of deep inelastic scattering at high energies is a popular 

and successful way of describing and analysing the HERA data on electron-proton 
scattering, that is in particular quasi-real and virtual photon-proton scattering. In 
the recent past it has turned out that the framework of the dipole picture is well 
suited for studying the region of relatively small photon virtualities Q^, that is the 
interesting transition region from perturbative to nonperturbative QCD. Particular 
interest concentrates on the problem of finding possible saturation effects occurring 
in that transition region. In essence, the dipole picture describes the deep inelastic 
photon-proton interaction as a two-step process in which the photon first splits into 
a quark-antiquark pair - a colour dipole - which then in the second step scatters off 
the proton. The first foundations on which the dipole picture rests were already laid 
in [5l [6]. For a general discussion of the dipole picture and for further references see 
[7]. One key problem in applying the dipole model in the transition region is to find 
a suitable description of the second step of the reaction, in other words to find the 
correct behaviour of the dipole-proton cross section. Models for this cross section have 
been developed and tested in inclusive as well as diffractive scattering processes, see 
for example [8]-[24]. Particularly popular is the Golec-Biernat-Wiisthoff model [8]-|10j 
which already in its first and simplest version gave a surprisingly good description of the 
HERA data. An important improvement of that model was obtained in [25] due to the 
inclusion of the correct logarithmic scaling violations according to DGLAP evolution 
at large while preserving the success of the model at low . More recently, also 
the problem of describing the dependence of the dipole cross section on the impact 
parameter and hence the question of the t-distribution of the cross section has been 
addressed [26l[27|. 

The two-step nature of the scattering process described above is most easily vi- 
sualised in a perturbative situation, and much of the theoretical development of the 
dipole model has been guided by the perturbative picture of deep inelastic scattering 
in QCD at high energies. Loosely speaking the dipole model tries to approach the 
transition to the nonperturbative region of small momenta from the perturbative side. 
Although its phenomenological success is rather encouraging it is not a priori clear that 
this method is in fact admissible. It is therefore very important to study whether there 
can be contributions to the cross section at low which are not compatible with the 
formalism of the dipole model, that is contributions which cannot be correctly accom- 
modated despite the freedom in modelling the low-Q^ dipole-proton cross section. It 
is one aim of the present paper and the companion paper [28] (hereafter referred to as 
II) to address this question. One of our main results will in fact be the identification 
of two contributions to the amplitude which are not contained in the simple dipole 
picture. Both of these correction terms are small at large but can become large 
at small . A second question which we consider to be important and which should 
be addressed in a nonperturbative framework is the following. It is frequently asserted 
that some sort of saturation must occur in real or virtual photon-proton scattering 
at high energies due to the Froissart-Martin-Lukaszuk bound on total cross sections. 
However, this bound applies only to hadron-hadron scattering and a priori is not valid 
for a current-induced reaction (see for instance [7J). It is usually argued that the dipole 
formula allows one to consider photon-proton scattering in essence as a hadron-proton 
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reaction. But if one really wants to substantiate this type of reasoning and maybe some 
day derive an analogue of the Froissart-Martin-Lukaszuk bound for scattering one 
certainly should have a solid foundation for something like the dipole picture. It is our 
opinion that in order to achieve this goal one should not rely on perturbation theory 
alone. 

In this work we start from a description of the real or virtual photon-proton scat- 
tering process in a generically nonperturbative framework developed by one of us in 
|291 [30] . This framework allows us to find a classification of different contributions to 
the cross section and to identify the approximations and assumptions that are neces- 
sary to arrive at the dipole model. We are able to identify two corrections to the dipole 
picture which can become large especially at low photon virtualities. In addition we 
find that our approach is also well suited to study two important problems which in 
our opinion have not yet been sufficiently appreciated so far. The first problem is that 
in QCD a mass-shell condition for quarks does not exist due to confinement. This is 
clearly not expected to be a serious problem in perturbative situations. However, if one 
wants to study high energy scattering in the nonperturbative domain this problem has 
to be taken very seriously. We will not be able to solve this deep problem in the present 
papers, but we hope to provide a framework in which at least the effects of assuming 
a mass-shell condition for quarks can be studied in more detail. The other important 
issue to which we would like to draw some attention concerns the electromagnetic gauge 
invariance in photon-proton scattering. It is well-known that gauge invariance in gen- 
eral imposes very strong restrictions on the amplitudes and on their factorisation. This 
has for example been discussed for the case of diffractive p production at high energies 
in [31j- 111 II we find that important corrections to the simple dipole model are in fact 
required in order to ensure gauge invariance. It turns out that only at asymptotically 
high energies the simple dipole formula is obtained from the nonperturbative frame- 
work as a gauge invariant expression. This immediately raises the question whether 
presently accessible energies are large enough for the dipole picture to give a gauge in- 
variant description in a sufficiently good approximation. Our framework should make 
it possible to address this important question in more detail. 

The two papers are organised as follows. In the first paper we describe the gen- 
eral nonperturbative framework of our considerations and derive a classification of the 
contributions to the process of real or virtual Compton scattering on a proton. After 
discussing the relative importance of the different contributions we proceed to study 
in detail those which we expect to be leading in the high energy limit. In the further 
discussion we then focus on these contributions and discuss in particular the renormal- 
isation of the photon-quark-antiquark vertex. The high energy limit of the amplitude 
obtained so far is then studied in the second paper. There we discuss how the usual 
formula for photon-proton scattering in the dipole picture emerges. Also in the second 
paper we discuss the important issue of gauge invariance. In particular we study how 
gauge invariance is manifest in the leading contribution to the scattering process. We 
then show that in the high energy limit different terms in the amplitude become gauge 
invariant separately. The perturbative wave function of the photon is derived and the 
consequences of gauge invariance for its correct definition are pointed out. Finally, we 
discuss some properties of the wave function of the photon and derive some simple 
phenomenological consequences of the dipole picture which can help to determine the 
limits of its applicability. We present our conclusions in II and point out some questions 
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that in our opinion deserve further study. Some technical details of our calculations 
are presented in the appendices of the two papers. 

2 Nonperturbative techniques and the dipole picture 

Here we analyse photon induced reactions using nonperturbative techniques introduced 
in |29l I30j . The aim is to see if one can justify the dipole picture for high energies or if 
modifications of the dipole picture are necessary. Our method is quite general and can 
be applied to many processes. For definiteness we consider real or virtual Compton 




pip') P{p) 
Figure 1: Compton scattering on a proton 



scattering on a proton (see figure [T]) : 

j(*)iq)+p{p) l^*\q')+p{p'), (1) 

where we indicate the four-momenta in brackets. For q' = q the absorptive part of the 
Compton amplitude can be related to spin-dependent and spin-independent structure 
functions. In the present papers we always work in leading order of the expansion 
in powers of the electromagnetic coupling aem- Note that throughout this paper the 
figures are drawn such that the incoming particles are on the right hand side. We 
have chosen this ordering in order to make the correspondence of the figures with the 
equations as close as possible - despite the fact that the opposite choice is probably 
more common in the literature. We suppose 

q^ = -Q^<0, q'^ = -Q'^<0, (2) 
but do not require = Q'^. The usual Mandelstam variables are 

t = (p-p'f = (q'-q)\ (3) 
u = {p- q'f = {p' - qf . 

Our procedure is rather general and allows us to treat the Compton amplitude for the 
case of real and virtual photons on the same footing. In the case of virtual photons 
our calculations apply to transversely as well as to longitudinally polarised photons. 
However, for the latter there are some subtleties in the high energy limit, as we will 
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explain in II. We will at first keep to the cartesian indices ;U, i/ = 0, 1, 2, 3 for the photon 
polarisations. 

Before going into the technical details let us now give a short non-technical descrip- 
tion of our method. We use the reduction formula to represent the amplitude for our 
process ([T|) as an integral over vacuum expectation values of the currents and the fields 
for the protons. Then these vacuum expectation values are expressed as functional 
integrals in the usual way. Since the Lagrangian of QCD is bilinear in the quark de- 
grees of freedom, q and we can perform the integrations over q and q. This allows 
us to give a classification of contributions to our amplitude in diagrammatic language 
according to the topology of the quark line skeleton, see figure [2] below. There the 
quark lines correspond to full quark propagators in a given external gluon field and the 
shaded blobs indicate the functional integration over all external gluon fields with a 
measure given below. The usual expectation is that at high energies the contribution 
shown in figure [2^ will dominate. We show that also the diagram of figure Wp can give 
a leading contribution at high energies. We study first figure [2^ in detail. We 'cut' 
the quark lines after the vertex where the photon 7(5, v) enters the diagram and insert 
suitable factors 1, given by the Dirac operator times the free Green's functions of the 
quarks. Integration by parts and the representation of the free Green's functions of the 
quarks as dyadic products of quark and antiquark wave functions lead to four terms 
one of which can be interpreted as expected in the dipole picture. The photon ^[q, v) 
splits into a quark- antiquark pair, which then interacts in all possible ways with the 
proton, see figures [TH and [8^ below. However, these quarks and antiquarks are still off 
their energy shell, and one has to integrate over all these ofF-shell energies. But at high 
energies the leading contribution comes from the on-shell quarks. This corresponds to 
the contribution of a certain pinch singularity in the amplitudes considered as functions 
in the (off-shell) energy plane. Adding now the contribution of the diagram of figure 
Wp, treated in a similar way, we have indeed described the original amplitude for the 
process ([TJ at high energies as a two step process. The leading term of the amplitude 
is given as a convolution of the bare ^qq vertex function with a scattering amplitude 
for a bare qq pair on the proton. It remains to introduce the renormalised quantities 
in place of the bare ones. This we do by using the Dyson- Schwinger equation for the 
^qq vertex function which contains all higher order QCD corrections to the vertex. 
Here our approach is inspired by the standard treatment of overlapping divergences in 
QED. Finally, we arrive at the expression shown graphically in figures [TO] and [TT] below. 
The leading contribution to the amplitude for the reaction ([1]) at high energies is given 
by two terms. The first one is a convolution of the renormalised ^qq vertex function, 
describing the splitting of the photon to the qq pair, with the renormalised amplitude 
for qq + p ^ iW + P- The second term describes rescattering corrections. Thus, 
on the side of the colour dipole amplitude, the contribution T^^^ from the diagram of 
figure [2)3 must be added, and on the side of the incoming photon, we find the rescat- 
tering correction. These additional terms are also leading at high energies and could 
be important for phenomenology especially at small photon virtualities. 

In addition to deriving this general structure of the amplitude for ([1]) at high ener- 
gies we discuss in II the problem of electromagnetic gauge invariance and the related 
subtleties for the amplitudes of longitudinally polarised virtual photons. In summary, 
from the contribution of figure [2^ and [2)3 we find something like the dipole picture at 
high energies. But usually only the term T^°) from the diagram of figure [2^ times the 
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photon vertex function is discussed in the dipole model. We find important additional 
terms as mentioned above. 

A key step in our procedure is to cut the amplitude next to the vertex of the in- 
coming photon and to insert a factor 1 consisting of the Dirac operator and of Green's 
functions of the quarks. A similar step has of course been used more or less explicitly 
in many studies of the Compton amplitude or of diffractive photon-proton scattering. 
Many of those investigations are, however, restricted to the case of virtual photons 
in which perturbation theory can be applied. In our method this step is used in the 
framework of a nonperturbative description of the real and the virtual Compton am- 
plitude. But several technical aspects of our method clearly resemble those studies, see 
for example [32l [33] and also [Ml ES] . 

After this outline of our procedure we now turn to the actual calculations. These 
involve sometimes lengthy formulae. But we hope that the figures accompanying them 
will, nevertheless, make them transparent. We should point out again that in order to 
make the correspondence between the equations and the figures as easy as possible all 
figures are drawn such that they are to be read from right to left. 

2.1 Classification of contributions to real or virtual Compton scatter- 
ing 

The matrix element for the reaction ([!]) is 

M>:r,{p',p,q) = ^ I d'^x e-'''^p{p',s')\T*J^^{0)r{xMp,s)) , (4) 

where nip is the proton mass, s,s' = ±1/2 are the spin indices, and T* is the covari- 
antised T-product. Our proton states are normalised to 

{pip', s'Mp, s)) = 2p\27Tf5,,s S^^Kp' - P) . (5) 

All our conventions regarding kinematics, Dirac spinors etc. follow |36j . 
The hadronic part of the electromagnetic current is denoted by 

j\x) = Y,Qix)Qil^Qi^), (6) 

where q = u,d, . . . are the quark field operators and Qu = 2/3, Qd = — 1/3, . . . are the 
quark charges in units of the proton charge. 

Let ijjpix) be a suitable interpolating field operator for the proton. We can take ipp 
to contain three quark fields, 

7pp{x) = raf3yUa{x)ui3{x)d^{x) , (7) 

and accordingly 

■0p(x) = taf3^d^{x)ui3{x)Ua{x) , (8) 

where Ta/3'y and Fq,^^ are coefficient matrices and a, (3, 7 summarize Dirac and colour 
indices. Explicit constructions of field operators ipp{x) of the type ([7]) can be found in 
[37 1 138 1 [39]. Let Zp be the proton's wave function renormalisation constant defined by 

{0\M^Mp,s)) = y%e-^P^Us{p). (9) 
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Using now the LSZ reduction formula [30] and the representation of Green's func- 
tions by the functional integral we get 



27rmpZp 



'y' +mp) 



(10) 



Here we use a shorthand notation for the functional integral. For any functional 
F{G,q,q) of gluon and quark fields we write 



{FiG,q,q)) 



G,q,q 



^ I ViG,q,q)F{G,q,q)exp 



d'^xCQCT){x) 



with 



Z = 'D{G,q,q)exp 



d'^xCqcDix) 



(11) 



(12) 



Of course, on the r.h.s. of pOj) . the quark-field operators in J^^ and tpp are to be replaced 
by the corresponding Grassmann integration variables of the functional integral. Gauge 
fixing and Faddeev-Popov terms in the functional integral are implied and not written 
out explicitly. 

In pop we have used translational invariance to shift the argument of the current 
to zero. In the following this makes our procedure look somewhat asymmetric with 
respect to the treatment of the two currents and J'^ . We could equally well make the 
treatment completely symmetric with respect to J'^ and J*^ by shifting the argument 
of ^pp or ^pp to zero in (fTO]) . see appendix El 

Since the Lagrangian of QCD is bilinear in the quark fields, we can immediately 
perform the q, q functional integral. This leads to a number of terms representing all 
possible contractions of quark and antiquark field operators, as is explained in ap- 
pendix El The resulting terms can be classified in terms of the quark line skeleton in 
the graphical representation of the amplitude. The different types of contributions to 
the amplitude are shown pictorially in figure [2l the explicit formulae for all terms are 
given explicitly in appendix E In figure [2] each quark line corresponds to a full quark 
propagator in a given external gluon potential. An integration over all gluon poten- 
tials with a functional integral measure including the fermion determinant is implied. 
Corresponding to the diagrams of figure [2] we get the following decomposition of A^^,*^ 

M^;^{p',p,q) = M^;l'"'{p',p,q) + . . .+Mi',l'"'{p',p,q) . (13) 
Putting all factors together we get for instance for ■M^^l^'^ (see figure [2ti) 



M'~fr{p',p,q) 



2TTmpZp 



dS dSe''"'uAp'){-i^y' +mp) 



(14) 



E Ql{ My')i^piy) ^^'^^^(9))^ i^}y +m,)^x,(p)e-^^ 



where the contraction tppiy')tpp{y) is defined in (lA.llh and 



d'^xTr 



7^4'^(0,:E;G)e-^«-7"5)?'(^,0;G) 



(15) 
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(a) (b) 




Figure 2: The diagrams resulting from the functional integral (lA.Sjl corresponding to 
, . . . , and to the amplitudes M^"^ M^^) pf 
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Here Sp\x,y;G) is the propagator for quark flavor q in the external gluon potential 
G, 

(z7'^D^-mW)4^)(a;,y;G) = -5(^)(x-y). (16) 

In ()14p we denote by {■)g the functional integral over the gluonic degrees of freedom, 
including the fermion determinant. That is, we define for any functional F[G] 



1 



G 



Z' 
exp 



- / V{G)F[G] JJdet \-i{n^Dx 



ml°) + ie] 



d'x^TT(Gxp{x)G^'{x)'^ 



where Z' is the normalisation factor obtained from the condition 

(1)g = 1. 



(17) 



(18) 



Again gauge fixing and Faddeev-Popov terms are implied. 

In a similar way we get from ()A.9P and ()A.12p for the amplitude ■M-^J^f''^ correspond- 
ing to figure 



-m2^(p',p,9) 



27n7ipZp 



(19) 



^Qg'Q,(V'p(y')V'p(2/)(-l)Tr 



7^'4''^(0,0;G) 



1 ,1 



d'^xTv 



YS^p\x,x;G) ) Ji^y +mp)us{p)e 



ipy 



The corresponding expressions for the other contributions A^*-*^-* - M.^^^ are given ex- 
plicitly in appendix lAl 

The methods used here for the classification of different contributions to the ampli- 
tude can also be applied to other current induced reactions at high energies. In the case 
of exclusive reactions for example one chooses appropriate interpolating field operators 
for the produced hadrons and can proceed exactly as discussed here for the Compton 
amplitude. 



2.2 Relative size of different contributions 

With the diagrams of figure [2] we have given a classification of all contributions to the 
Compton amplitude, see (I13p . At high energies ^/s we expect the diagrams of figures 
[2ji and [lb, that is M^"^ of ^ and M^^^ of to give the dominant contribution. 
Indeed, let us consider the diagrams of figures [2^1-12^ for large ^/s. In figure [2t for 
instance we find at Born level that a large momentum of order ^/s has to fiow through 
the upper quark line. This leads to a suppression of order l/\/i relative to the diagrams 
of figures [2^, [2b. Similar arguments apply to the diagrams of figures [2ll-[2^. 

Going one step further in sophistication we can have a look at resummed per- 
turbation theory. In the perturbative framework the generic diagrams of figures [2^ 
and [2b include gluon ladder diagrams, see figures [3] and [U respectively. These are, 
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Figure 3: Perturbative gluon ladder diagram which is contained in the generic diagram 
of figure for the amphtude Ad^""^ 




Figure 4: Perturbative gluon ladder diagram which is contained in the generic diagram 
of figure [2b for the amplitude M ^''^ 



of course, just the diagrams containing the leading order BFKL pomeron exchange 
|4H H2] . The leading perturbative diagrams contained in the generic diagrams of fig- 
ures[2}:-[2^, on the other hand, do not contain such gluon ladders. Instead, they typically 
contain quark-antiquark exchanges in the t-channel. As an example figure [5] shows a 
perturbative quark ladder diagram contributing to the generic diagram of figure [2}:. 
But it is well known that diagrams with quark-antiquark exchanges (or more generally 
fermion-antifermion exchanges) in the i-channel are suppressed by powers of the energy 
with respect to the corresponding diagrams containing gluon (or more generally vector 
boson) exchanges. Similarly, also the generic diagrams in figures [2]1 and [2^ contain 
quark-antiquark exchanges as the leading terms in perturbation theory. Presumably 
also the leading terms of the generic diagrams of figures Ef and [2^ correspond to quark- 
antiquark exchanges, see figure [H where this is shown for the diagram of figure [2^. We 
can hence conclude that the leading perturbative diagrams for the generic diagrams of 
figures [2]3-[2^ are suppressed by powers of ^/s with respect to the leading perturbative 
contributions to the amplitudes Ai^""^ and Ai^^\ 
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I_l '\/\/\/\y\y\j^ - 




Figure 5: Perturbative quark ladder diagram which is contained in the generic diagram 
of figure [2]:: for the amphtude Ai^'^^ 




Figure 6: Perturbative diagram which is contained in the generic diagram of figure [2^ 
for the amphtude TW^^^ 



We also arrive at the same conclusion when we switch from perturbation theory to 
Regge theory. Quite generally, we expect the dominant contribution from the diagrams 
of figures [2^1 and [2)3 to be due to pomeron exchange, soft and hard, with its well-known 
phenomenological properties, see for instance [7]. On the other hand, we estimate that 
the main contribution from the generic diagram of figure [2t for large ^/s is due to 
reggeon exchange which would correspond to the diagram of figure [5l now interpreted 
as a diagram of Regge theory. Again, we know from phenomenology that reggeon 
exchanges are suppressed by a factor of roughly l/^/s relative to soft pomeron exchange, 
see [7\ . Thus we can argue that the generic diagrams of figure [2}: and in a similar way 
those of figures EJi and [2^ correspond to reggeon exchanges and are thus not leading 
for large ^/s. The same holds for the diagrams of figures [2]F and [2^, as is illustrated in 
figure m Again we find that the amplitudes Ai^'^^ - A4 are subleading at high energy 

VS. 

Let us now discuss the relative magnitudes of the amplitudes A^^*^-* and Ai^^'^ cor- 
responding to the figures [3] and [H At high energies the leading contributions to both 
amplitudes have the same dependence on the energy. As discussed above that depen- 
dence results from Pomeron exchange or from gluon ladder diagrams, depending on 
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whether one describes it in terms of Regge theory or perturbative QCD (pQCD). The 
typical perturbative diagrams are shown in figures [3] and HI respectively, and in pQCD 
these diagrams represent the leading contribution at high energies. The main difference 
between the two amplitudes is the coupling of the gluon ladder (or the Pomeron) to 
the incoming and outgoing photon, or in perturbative language the impact factor. In 
the impact factor in A^(*) there are two quark loops which are connected by gluon 
lines. Due to the charge parity of the photon a minimum number of three gluons is 
required to couple to each of these loops, and the diagram shown in figure H] is hence the 
lowest order contribution to the impact factor of Ai . Compared with the amplitude 
A^(") the amplitude Ai^'^^ therefore contains two additional factors of ag which are not 
compensated by any large logarithms. The relevant scale of ag is determined by the 
photon virtuality Q^. Consequently, the amplitude Ai^^^ is strongly suppressed at high 
virtualities. At small virtualities, however, Og becomes large and that suppression is 
no longer relevant. Besides these factors of there appears to be no other simple 
mechanism suppressing the diagrams of figure HI Accordingly, these diagrams can po- 
tentially give an important contribution in the region of low photon virtualities. In our 
opinion their contribution, which is often neglected in applications of the dipole picture 
of high energy scattering, deserves further study. It would be especially interesting to 
find ways of estimating the relative size of the terms M^"^ and at low 

usmg 

experimental data. 

It should be noted here that in our development of the dipole picture below we 
find that perturbative diagrams contained in figures [2ti and [2)3 contribute to the colour 
dipole amplitude side and photon wave function side such that one can speak of a 
mixing of the diagram classes (a) and (b) in the colour dipole picture. This is relevant 
in higher orders in Og as we will discuss when we make the transition to renormalised 
quantities in section 12.41 below. The preceding discussion refers to the leading terms 
(at large ^/s and large Q^) and is therefore not affected by this effect. 

From the observations made here we conclude that at high energies the amplitude 
Ai^'^^ is expected to give the leading contribution if the photon virtuality is large, and 
that at lower virtualities also the amplitude Ai^^^ should become important. We will 
therefore in the following discuss in detail the amplitude A4^"'\ given in (|14p . which 
corresponds to the diagram of figure [2^. For the amplitude Ai^^\ figure [2)3, the same 
steps can be done in a completely analogous way, see appendix ICl 

Finally we note that for the hypothetical case of scattering of an isovector photon 
on an fi" target only the diagram of figure [2^ exists as long as one assumes equal 
up and down quark masses and neglects higher orders in the electromagnetic coupling 
Uem, as we do throughout the present papers. Indeed, in this case the current has the 
structure 

JL.ector = 1^7'^ -^d^d, (20) 

whereas the consists of three s-quarks. Thus none of the diagrams of figures 
can be drawn for this case and A^^"^ gives the complete amplitude. 

2.3 From current to quark scattering for a fixed gluon potential 

In the following we consider the diagram of figure [2^1. We will in essence cut the quark 
lines at the vertex marked x, insert factors of 1, and relate the amplitude for the 



11 



scattering of the photon on the proton to that of scattering of a ggpair on the proton. 
We recall first that the free quark propagator for mass nig satisfieq^ 



-- 6^^\x-y). 



ly +mq) 



(21) 
(22) 



A simple exercise shows that S"^'^^ can be written as a sum of dyadic products of quark 
and antiquark wave functions, see (jB.lOp of appendix [Bl 

We shall now use ^ and ([BJO]) to rewrite A(«)^''(^) of ([ED in a form where 
we can make the transition from current scattering to qq scattering. Inserting factors 
of one in (fT5|) we get 

Ai<i)t''^{q) = j (fx (fz d^z' Tr [t^^^^^O, z; G)5^*'^ {z - x) 
e-'^^^^S^^\x- z')S^^\z',Q-G) 



(fx (fz (fz' Tr 



^^^Sf{^,z-G){-i^,+m,)Sp 
■*'?^7"4''°^(^> ^'){^z' +m,)sf{z\ 0; G) 



(9,0) 



Z.X] 



(23) 



Integration by parts and use of (jB.lOp leads us to a representation of A^''^^'^{q) as a 
sum of four terms which we can interpret in terms of q and q scattering. The key 
step in using (|B.10p is a decomposition of the free propagator S^p'^^ as a spin and 
colour sum over Dirac spinors, see (|B.4p . (jB.Sp . Those spinors immediately lead to an 
interpretation in terms of quarks and antiquarks. The details are given in appendix iBl 
The final result derived there is (see ()B.25P ) 



where 



with 



ai 



2-K J Lo + ie J (27r)32A;0 



■aj, 



-^(fc5'),r'|i?('?)'^|A:^,r)n,(A:)7^r;,,(fc'«) 

r,r' 

(gO_^o_^/{i)o_^ + .g)-i(2fc/(i)0)-i^ 



(24) 



(25) 



(26) 



«2 = j;(A;:(2)^^/|^{g)M|^^^^)^^,(^)^.^^,(^/{2)) 
r,r' 

(_^o ^ ^0 _ ^/(2)o ^ ^ ^ ie)-i(2A;'(2)0)-i , (27) 

^Note that the superscript in the quark propagator S^'"' indicates that this is the free propagator, 
whereas in the quark mass m'^^ it indicates that this is the bare mass, see (|A.1|) . 
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«3 = -5^(A:5=^),r'|i?('?)'^|A;^,r)t;,(A:)7^n,,(fc'(3)) 

r,r' 

(_gO _ ^0 _ ^/(3)0 _ ^ ^ ie)-i(2A;'(3)0)-i , (28) 



a4 = 5](A:5'\r'|i?(5)^|fc^,r)t;,(A:)7'^^;,,(fc'(')) 

r,r' 

(gO + ^0 _ ^/(4)o ^ ^ ^ ie)-i(2fc'W0)-i . (29) 

The operator 5^'')'' is defined in ([Bl4]) and k'^^^k'^^ in ([R26l) - ([B:29l) . Here r,r' are 
spin and colour indices with r = (A, A), X = ±1/2 and A = 1, 2, 3. The four terms are 
illustrated in figure [71 In that figure the dashed vertical line indicates an integration 
over uj and over the three-momentum k according to (|25p . We emphasise that the 
vertical line does not imply that all quark and antiquark lines are set on-shell, see 
the discussion below. The arrows on the lines indicate the interpretation as quark or 
antiquark. The orientation of the momentum assignment to the lines is from right to 
left for all lines, irrespective of the arrows on the lines. The interpretation of A^^^^'^ (q) 
is now easily obtained (see figure (T]) : 

• The term A^^^^'^{q) corresponds in essence to the splitting of the initial real or 
virtual photon into a qq pair with momenta k,k'^^\ This is combined with the 
amplitude for a qq pair of momenta k^, k'S^^ scattering on the fixed gluon potential 
and recombining to the final real or virtual photon q' (figure [7h) . 

• The term A^2^^'^{q) corresponds to a quark of momentum /c'^^^ absorbing the 
photon q and going to a quark of momentum k. This is combined with the 
amplitude for a quark k^j shaking off a photon q' in the fixed gluon potential and 

/(2) 

going to quark k^j (figure Ub) ■ 

• The term A^^^^'^lq) corresponds to an antiquark k and quark k'^^^ annihilating by 
absorption of a photon q. This is combined with the amplitude for the creation 
of a quark k^j , antiquark k^^ and photon q' in the fixed gluon potential (figure 

Eb). 

• The term A^^^'^ (q) corresponds to an antiquark k absorbing a photon q and going 

to an antiquark k'^'^\ This is combined with the amplitude for an antiquark k'S^^ 
going to an antiquark fc^ with emission of a photon q' in the fixed gluon potential 
(figure EH). 

Several points are noteworthy here. First, we have nowhere made the assump- 
tion that free quarks exist. Second, the mass parameter niq in the above formulae 
is completely arbitrary. Further, the spinors u and v in the splitting of the pho- 
ton q into a quark-antiquark pair etc. and the spinors in the scattering amplitudes 

{k'S^\r'\B^'^^^\kui,r) etc. are always those of quarks on the mass shell niq, see (jB.Sp . 
But in the plane wave factors on the r.h.s. of ()B.8P the energies 5 kijj etc. are not 
equal to the energies of on-shell quarks, see ()B.7p and ()B.26P - (]B.29p . For the diagrams 
in figure [7] this means that the quark and antiquark lines to the right of the vertical 
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Figure 7: Interpretation of the four terms Aj (q) of ()25p where j = 1, ... ,4 corre- 
sponds to (a), . . . , (d). The shading indicates the fixed gluon potential. The arrows 
on the hues indicate the interpretation as quark or antiquark. The orientation of the 
momentum assignment to the Hnes is from right to left for all lines, irrespective of the 
arrows on the lines. 
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line have energies corresponding to on-shell quarks. The quark and antiquark hnes 
to the left of the vertical line, on the other hand, have energies corresponding to off- 
shell quarks, and these energies contain the integration parameter to. In this respect 
our result is reminiscent of noncovariant perturbation theory. But we stress that no 
perturbation expansion has been made in our approach. 

We should add a remark here concerning the interpretation of the different terms 
in figure [71 On first sight the occurrence of the last three diagrams might appear 
counterintuitive. We emphasise that all diagrams are obtained from the term 
corresponding to figure [2^. The decomposition into the four terms of figure [7] originates 
from the decomposition of propagators into spin sums. In that sense the free propagator 
contains all possible quark and antiquark spinor states, giving rise to the different 
combinations of incoming and outgoing quarks and antiquarks in the four terms in 
figure [71 Also that decomposition reminds one of noncovariant perturbation theory. 
Loosely speaking, the quark loop of the upper part of figure [2^1 is in figure [71 cut and 
deformed in such a way that the intermediate quark and antiquark are 'bent' into the 
initial or final state. 



2.4 From current to quark scattering in QFT 

Now we are ready to insert our result for ((Ml)- ([29]) in 1^ and to analyse the 

resulting expression. We get 

M'^l^''ip',p,q) = Y.M^;f^''ip',p,q) , (30) 



Mi':f'"'ip',p,q) = -^£-^Jdyd^ye'P'y'uAp)H^y'+rnp) (31) 

9 



uj + iej (27r)32A;0 
(2A;'«0)-ij;(7(9',/x),p(p',s')|T('^)k-(A;L(i\rO,(?(fe.,r),p(p,s)) 

Ur{k)ZgYVr'{k'^^^) , (32) 



uj + ieJ (27r)32A;0 
{2k'^'''>^)-^Y.^^{q',^,),q{k'^^\r'),p{p\s')\r^-^^^^^ 

Ur{k)ZqYUr'{k'^^^) , (33) 
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(2t'IW)-i^{7(<,',M),f(fc„,r),5(t-i",r'),p(p',s')|rl»>b(p,s)> 

r',r 

i;,(A:)Z,7V'(^'^'^), (34) 
(2t'l''l»)-i^(7(<(',M),9(*^,r),p(p',/)|rl°l|f(til''l,r'),p(p,s)> 

{-l)Vr{k)Z,^^Vr'{k'^''^), (35) 

where we have defined 

(7(9',;"),p(y,5')l^^''^k"(A:5'\^'),'?(^-,0,P(P,s))= (36) 



V'p(y')?p(2/)(-l)(fc^^'\^'|S^'^1A;^,0)^ (^^, +"ip)^.(p)e-'^^ 
(7(9',/^),9(A;5'\rO,M/,^')l'^^"^k(fc-,0,P(P,s))= (37) 



ll:jXy')i^^{y){k'P,r'\B^'i^^^\k^,r)) {t^ ^ +mp)n,(p)e-*^'^ 

G 



= - W / ^^'y' Ae^^'^'n,,(y)(-^?y +mp) 



S^^piy)i-l)ik'S'\r'\B^'i^^\k^,r)) (z^y +mp)u,{p)e-'P\ 



M,^^),qiku.,r),pip',s')\T(''^\qik'S''\r'),p{p,s))= (39) 

dydSe'P'y'uAp'){-i^y, +mp) 



27r?71p Zq 



^^V'p(y)(-l)(A:5'\r'|i?(^)'^|fc^,r)) (i^, +mp)n,(p)e-^P^ . 



Here it seems that we can already make a connection to the usual dipole picture for 
photon-induced reactions. The amplitude (j32p has been written as a product 

of a term describing the splitting of the initial photon ^{q, v) into quark and antiquark 
and the subsequent scattering of the quark and antiquark off the proton to give the 
final state \p{p' ■,s'),^{q\^)). But there are two obvious complications. 

First, we have written }A^°'^ as a sum of four terms in (j30p . The interpretation 
of these terms is easily obtained from ([32|) - ([35|) and is shown in figures [8^-[8ll. It 
is analogous to the interpretation of the four diagrams of figure [7| in the preceding 
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(a) 




Figure 8: Interpretation of the four terms M^^'J\p' ,p, q) ([32|) - ([35|) where j = 1, . . . ,4 
corresponds to (a),. . . ,(d). The shading indicates all possible interactions of the proton 
with the quarks at the top of the diagrams. The arrows on the lines indicate the 
interpretation as particle or antiparticle. The orientation of the momentum assignment 
to the lines is from right to left for all lines, irrespective of the arrows on the lines. 
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section. The vertex diagrams on the r.h.s. of figure [8] correspond to the vertex factors 
Ur{k)Zqj'^Vr'{k'^^^) ctc. in ([32|) - ([35|) . The diagrams on the l.h.s. of figure [8] correspond 
to the matrix elements (|36p -(|39 p . We emphasise that these are matrix elements for 
quarks off the energy shell. The normalisations in (|36 p - (|39p are, however, chosen with 
correct LSZ factors and wave function renormalisation constants in such a way that the 
matrix elements of T*^") are finite and would correspond to physical T-matrix elements 
under the following conditions: 

• The quarks q are supposed to have a mass shell with a pole mass mg^poie and we 
choose our mass parameter - which is still at our free disposal - equal to the pole 
mass. That is, we set nig = mg^poie- 

• We replace the momenta of the quarks and antiquarks in (I36|) - (l39p - which are 
in general off the energy shell - by the on-shell momenta. That is we make the 
replacement k^^ —>■ k, k^^ —>■ k^^^ etc. 

Keeping these points in mind we see that figure [8^ is clearly interpretable in the dipole 
picture, in contrast to the other three diagrams. Figure [SJa describes the absorption of 
the photon g on a quark k'^'^^ going to quark k. This is combined with the amplitude 
for scattering of quark fc^ on the proton giving quark k^ , a photon q' and the proton 
p' . The interpretation of the figures Et and [Hi is analogous. 

A second complication arises when we consider the renormalisation of the photon 
vertex. We have inserted a wave function renormalisation factor in the scattering 
amplitudes T^"^ (I36|) -(i39]). since otherwise they would not be finite. Correspondingly an 
explicit factor Zg has been put into the splitting function of the initial photon ^(q, u). 
Thus we are still dealing with (in general) divergent integrands in (j32l )-(l35]) due to the 
explicit factor Zg. This divergence in Zg must, however, be cancelled by a divergence 
in the integral over uj and k. 

The situation here is similar to the case of overlapping divergences in QED. The 
way to deal with such overlapping divergences is well known (see for instance |43j ) . 
Let r('^)'^(/, be the renormalised ^qq vertex function. This function satisfies a Dyson- 
Schwinger equation shown diagrammatically in figure [9l which is obtained as a straight- 
forward generalisation of the corresponding equation in QED, see i3]. Here JC^'''''') is 
the renormalised kernel of q'q' to qq scattering which contains all connected Feynman 
diagrams which are quark-antiquark two particle irreducible in the s-channel and S^p ^ 




Figure 9: The equation satisfied by the renormalised ^qq vertex function in diagram- 
matic form. In the last term a sum over the quark fiavour q' (see eq. (|40p ) is implied. 
The arrows on the lines indicate the interpretation as quark or antiquark. 
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is the renormalised quark propagator. Since we neglect all higher order terms in the 
electromagnetic coupling Oem the diagram expansion of the kernel comprises only QCD 
diagrams. Among the diagrams contained in is the one-gluon exchange between 

the quark and antiquark, as well as the annihilation of the quark and antiquark into 
three gluons in the s-channel which then create another quark-antiquark pair (compare 
the diagram in figure [T3l below). In the latter type of diagram the photon can couple 
to a qq pair of a different flavour than the one on the left. This possibility requires the 
two flavour indices in the notation K^'^''> ^ . 

In a shorthand notation the Dyson-Schwinger equation reads 

r('?)^(/, /') = Z^Y K'^'^'"''^ S^fT'^'i> S^f . (40) 

q' 

Solving for ^^7^ we get 



(41) 



All quantities on the r.h.s. of this equation are renormalised quantities, in particular, 
they are finite. The divergence of Zq originates from the integration in the last term. 
Replacing now Zq^'^ in (j32p - (l35l ) according to (I4ip we get expressions for A^^"'-') for 
j = 1, . . . , 4 where all integrands are finite. We still have some freedom how to choose 
the momenta V on the r.h.s. of ()4ip . A natural choice for is 

l = k, I' = -A;'(^) . (42) 

With this the final result for A^^"'^-* is shown in diagrammatic form in figure fTOl and 
reads 

(2A:'«0)-i^(7(<?',/i),p(p',/)|r('^)|g-(A.5i),r'),g(A;^,r),p(p,s)) 



(43) 



The explicit expressions and diagrammatic forms for A4^"''^\j = 2,3,4 are analogous 
and easily constructed. In II we shall study the high energy limit, where we will find 
that TW^"'^^ gives the leading term. 

The integrations over oj and k ()43p should be convergent and finite, but only for 
both terms in the curly brackets inserted together. If we insert for instance only r'-'^^'^ 
a divergence could occur. It would be very interesting to study in detail how the 
divergences cancel in this equation. It is likely that already a one-loop calculation 
would reveal the structure of this cancellation, possibly even in a simple model like the 
abelian gluon model proposed in [44j . 
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Figure 10: The two terms for the amplitude (j43|) in diagrammatic form. Again 

a sum over the quark flavours is imphed in the second term. The arrows on the lines 
indicate the interpretation as particle or antiparticle. The orientation of the momentum 
assignment to the hues is from right to left, irrespective of the arrows on the hues. 



For the amplitude Ai^^'^ corresponding to the generic diagram of figure [5)3 ah the 
same steps can be performed as for M^"''^. We get again a decomposition as in (j30p . 



(44) 



where of these four terms the important term at high energies is 

g",g 



1' J 



Vr' {k 



(45) 



Graphically this is represented in figure [TTl The explicit expressions are given in ap- 
pendix[Cj The matrix element {'y,p\T^'i ''^'^Iq^Q^p) in (|45p represents the scattering of 
a qq dipole on a proton giving a photon component q"^^q" and a proton. 



20 



Figure 11: The two terms for the amphtude (j45|) in diagrammatic form. The 

interpretation of the hnes is as in figure [TOl 



It is instructive to ihustrate the overlapping divergences problem which is solved 
by (j43p and (j45p with two simple perturbative diagrams. Consider first the diagram of 
figure [T2I Cutting this diagram along the dashed line gives a higher order correction 
to the photon vertex function. Cutting it along the dotted line we find a higher order 
correction to the scattering amplitude of the colour dipole on the proton. In the expres- 
sion (j43p both corrections are taken into account and the 'double counting' problem is 
exactly taken care of by the kernel term K^'^''^ \ 

This double counting or overlapping divergences problem in the 'cut' diagrams and 
its solution in (j43p and (j45p also 'mixes' the diagrams of figures [2^ and [2)3. This is 
despite the fact that these diagram classes are clearly identifiable and do not mix at the 
starting level. Indeed, consider figure [13l which is clearly a diagram of type [2)3, and cut 
the diagram first along the dotted line. Then we have a colour dipole-proton amplitude 
of type [2)3 with a simple vertex on the photon side. Cutting along the dashed line, we 
can interpret it as a colour dipole-proton amplitude of type [2^ with a vertex correction 
on the photon side. This shows that great care has to be exercised when cutting 
perturbative diagrams and trying to interpret them in terms of colour dipole amplitude 
times photon wave function. We would like to stress that the original classification of 
diagram types in figure [2] is not affected by this. The 'mixing' of diagram classes comes 
into play only when one considers cut diagrams, and then an individual diagram like 
the one in figure [13] can be associated with higher order corrections occurring to the 
left or to the right from the cut, depending on the position of the cut. 

Finally we want to make some remarks concerning the outgoing photon. As we 
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Figure 12: Perturbative diagram contributing to the amphtude M^"''^^ 




Figure 13: Perturbative diagram contributing to the amplitude but also to 



indicate briefly in appendix|A]one can treat the outgoing photon 7(g',/x) of the Comp- 
ton amplitude in complete analogy to the incoming photon 7(g, z^), performing all the 
steps we have done for the incoming photon in the same manner, including the renor- 
malisation of the the photon vertex via its Dyson- Schwinger equation. The resulting 
expressions would in their form be completely symmetric in the incoming and outgoing 
photon. For the purpose of the present papers, especially for the study of the high en- 
ergy limit to be discussed in II, it will however be sufficient to treat only the incoming 
photon in this way. 

3 Summary 

In the present paper we have studied the real and virtual Compton amplitude for 
photon-nucleon scattering in a nonperturbative framework. We have classified the 
different contributions to the amplitude according to their quark line skeleton. We 
have discussed the relative magnitude of the different contributions and have identified 
the two contributions A4 ^""^ and A4 represented by figures [2^ and [2]3 as the leading 
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ones at high energies. Concentrating on these two contributions we have then inserted 
on the quark hues after the incoming photon vertex suitable factors of one, given by the 
free Dirac operator acting on the free quark propagator for a chosen mass niq. After 
integrating by parts we have used spin times colour sum decompositions for the free 
quark propagators to obtain expressions which bring us closer to the dipole picture of 
high energy scattering. For both amplitudes A^(") and M^^^ we obtained four terms 
where in each case only one had the structure of a vertex function describing the 
splitting of the photon into a qq pair convoluted with a scattering amplitude for an 
off energy shell qq pair on the proton. The results for these parts of the amplitudes 
and M^^'^ are given in (jl3]) and ([M]), respectively, see figures [10] and [TTJ In II 
these results will allow us to isolate the leading term in the high energy limit. A key 
point in the present paper is that we have taken into account the renormalisation of the 
photon-quark-antiquark vertex. We find that the renormalisation induces an important 
rescattering correction which is usually not taken into account in the dipole picture. 

From the contribution of figure [2^ we will be able to obtain the usual dipole picture 
of photon-proton scattering after suitable further approximations. We will in fact 
show that this is the leading term at high energies and large photon virtualities. We 
find, however, that there are two important correction terms which can potentially 
become large at small photon virtualities, that is in the interesting transition region 
from the perturbative to the nonperturbative regime. The first correction comes from 
the rescattering of the quark and antiquark into which the virtual photon splits, the 
second correction comes from the contribution of figure [2)3. Both of them are suppressed 
only at large photon virtualities but not due to large energies. 

The motivation for the dipole picture of high energy scattering is closely related 
to the notion of light-cone wave functions, see for example I46j . When applied to 
a real or virtual photon this concept naturally suggests to include in addition to the 
quark-antiquark component also higher Fock states of the photon like for example the 
quark-antiquark-gluon component of the wave function etc., and to sum the amplitudes 
for the scattering of these components off the proton. Also the perturbative framework 
can be arranged according to this view. We emphasise that in our treatment nothing 
is missing and, in particular, higher Fock states do not and should not occur explicitly 
here. Recall in this context that the solution of the overlapping divergences problem 
for the vacuum polarisation in QED (see [43] ) also does not involve analogues of higher 
Fock states for the photon. Nevertheless we think that it is an interesting problem 
for further investigations to try to see where and how the various terms of a Fock 
state expansion of the photon as discussed for instance in [33] can be identified in our 
approach. 

In the second paper we will discuss in detail the high energy limit of the expressions 
(j43p and (j45p obtained here from the contributions of figures [2^ and [2)3. We will also 
study the problem of gauge invariance. In particular, we will discuss whether different 
parts of the amplitude are separately gauge invariant. We point out the consequences 
of this issue for the correct definition of the perturbative wave function of the photon. 
Finally we derive some phenomenological consequences of the simple dipole picture 
which can help to find its range of validity. With this brief outlook we want to close 
the present paper, further conclusions will be presented in the second paper. 
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A Propagators and decomposition of the Compton am- 
plitude 

In this appendix we give the derivation of the analytic expressions for the diagrams of 
figure [21 We start with the Lagrangian of QCD 

Cqcd{x) =-^Tv [Gxp{x)G^P{x)) + ^ q{x){i^^Dx - mf^)q{x) , (A.l) 

where G\p is the matrix-valued gluon field strength tensor and D\ the covariant deriva- 
tive (for the notation see |36j): 

Gap(x) = (^dxG;{x)-d,Gl{x)-g^'^UcGi{x)G';,{x))^, (A.2) 
Dxq{x) = (^dx + ig^'^^Gl^'^qix). (A.3) 

Here g^^^ is the bare coupling parameter and rriq^^ are the bare quark masses. Since 
^QCD is bilinear in the quark fields we can immediately perform the q and q functional 
integrals in (jlOp . To write the result in a concise form we use the quark propagator 
Sp\x,y;G) for given gluon potential, see ([TBI) . This propagator satisfies Lippmann- 

Schwinger relations. Let Sp''^^ be the free propagator for mass niq (see ([2T]) . (p2]) ). 
Then we have in matrix notation 

= s'-f^ - s'-^\g) (<7(°) r y + - m,) s}?'") . (A.4) 

From (1A.4P we get 

4^)(x,y;G)(f:^,+mW) = s'^^\x,y;G) (^i}^ +g(^^ + 

= 6^^Hx-y). (A.5) 
We define the contraction of two quark fields of flavour q as 

q{x)7!{y) = -S^^\x,y;G) (A.6) 
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and represent it by a line in our diagrams: 



X y 

In a functional integral over gluon and quark degrees of freedom we then have to 
contract all quark fields in all possible ways and integrate over the remaining gluon 
degrees of freedom including the fermion determinant. Schematically we have: 

(...q{x)q{y)...) = ( . . .-S^\x,y;G) . . .) 

\ I G,q,q \ I I G 

+ all other contractions of quark fields. (A. 7) 

Here the functional integral on the l.h.s. is defined in (jlip . ()12p . that on the r.h.s. in 
(|17p . (jlSp . For the functional integral in (jlOp (there written already for x' = 0) we find 
in this way 



i^,{y')J'{x')r{x)i^p{y) 

G,q,q 



= ^ Qq'Qq{^a'p''y'Ua'{y')up'{y')d^'{y')q {x')-i^q'{x')q{x)-i''q{x) 
q',q 

fa(i^d^{y)up{y)ua{y)) ^ 

I G,q,q 

= J^"^ + . . . + J^^^ . (A.8) 

The integration over quark and antiquark fields gives rise to a number of terms in 
which these fields are contracted in all possible ways. The resulting propagators can 
be understood as defining a quark line skeleton for the graphical representation of 
the respective term. We have classified the different types of diagrams in figure [2], 
and the terms J^""^ to J^^^ in the last line above are defined as corresponding to this 
classification. 

Inserting the expression (jA.Sh in ([TO]) leads to the decomposition of the amplitude 
in ([T3D with 

[ d^x e-^«" J^''^ (i ^„ +mp)us {p)e-'Py , (A.9) 

J x'=0 

and M^^^ to M^^) are related to J^''\ . . . , 

in a completely analogous way. The 
expressions x' = are represented in diagrammatic language in figure 

[21 The general rule is that for each quark line in a diagram in figure [2] one has to insert 
a factor —iSp\x, y; G) according to ()A.6p . The vertices with the photon just represent 
the coupling 'f'^Qq- The functional integral ( )g as defined in (fT7|) has to be performed 
over the integrand resulting in this way. This is indicated by the shading in figures [2^ 
to [2^. The explicit expressions are 



r\s^^\x',x;G)r-S'^\x,x';G) 



) , (A.IO) 

G 
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where 



{741'/3(y'>y;G)^4a'a(2/''y;G) -(«-/?)}, (a.h) 



further 



/ I 1 



q ,1 



TV 



I G 

(A.12) 



and 



-S^^^ {y', x'- 0)^1^ -Sf {x', X- G)Y-Sf {x, y; G) 
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]sP,,,iy',y-,G)]s^;l,^{y',y;G) 
+ Ql\sP^,,{y',y;G) 

-SP iy', x'- G)r-SP {x', x; G)Y-SP (x, y; G) 



- (a' ^ /?') - (a ^ /3) + (a ^ /3,a' ^ /3')|^ 
IS as but with 11^1/ and X ^ x'. Exphcitly, 



(A.13) 



G 



-S? iy', x; 0)7^^-4"^ (x, x'; 0)7^^-5? (x', y; G) 



lS';>M,y-,G)^S^^l^{y',y;G) - (a - /3) 



1 ^(n) 

+ Q«^41'^(2/'>y;G) 

i 



Sf{y\ X- Gyr-Sf{x, x'- G)r-SP{x\ y- G) 



4a'a(2/'-y;G) 



(a ^ /?') - (a ^ /3) + (a ^ /3, a' ^ /?') 



(A.14) 



G 
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We further have 



(-S'^;\y', X- G)Y-S'^\x, y- G)) (-Sf{y', x'- G)^^- 
- (a' ^ 0) - (a ^ /3) + (a' ^ /?', a ^ /?) 



i5?)(a;',y;G)') 



+ QuQa 



A * ^ / 7'7 

f x; 0)7^^-4") (x, y; G)] y; G) 

- {a ^ 13') - (a ^ /3) + (a ^ 0,a^ /?) 

V * ^ / 7'7 

(-SP{y', x'- G)AsP {x', y; G)) -SP^,^{y', y; G) 



(A.15) 



Y^Qg,i-l)TT (^j'^-^sf{x',x';G)^ 
\Qu-S^Fl,^iy', y; G) \ {-S^^\y', x; G)Y-SP (x, y; G)) -S^^^^iy', y; G) 

V L \ / a a 

- (a' ^ /?') - (a ^ /?) + (a' ^ a ^ /?) 
+ Qd(-S?{y', x; Gh''-S^^\x, y; G)) 

^ ^ , (A.16) 



G 
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{Q.-4y,(y'> G) \ (-SPiy', x'; Gh^-SP{x', y; G)] -S^^l^{y', y; G) 
- {a ^ 15') -{a^ 13) + {a ^ /?', a ^ /3) 
+ Qd i-Sf {y', x'- G)^^^-Sf {x', y; G)) 

i4"iJy^y;G)i4"j,^(y^y;G)-(a-/3)]|\ . (A.17) 

Finally, we indicate how we can make an analysis of the amplitude (|1U|) treating 
incoming and outgoing photons in a symmetric way. Using translational invariance of 
the functional integral we can write (jlOp also as follows 



2TTmpZp 



<fy(fx'(fx 



My'V Jn^')e"'"J''(^)My))^ A'^y +mp)us{p)e-'Py 

' G,q,q 



.(A.18) 



Inserting here the decomposition (jA.SP of the functional integral we obtain 



iTTnipZp 



d'^y d^x' d^x g^'j'^'e-^'^^e-^P^ 



Us'{p'){-i^y' +mp) J^''\i^y +mp)us{p) 
and similarly for M^''\ . . . , M^^) , With (|AlOl) we now find 



(A.19) 



y'=o 



M^fr{p',p,q) 



2'nmpZp 



I 2 



-ipy 



(A.20) 



j/'=0 



where 



d'^x d'^x Tr 



eii'^'^i'S^^\x',x-G)e-'"'^YSf{x,x';G) .(A.21) 



y 0(1)1 



The further analysis of A^''^^'^ can then be done exactly as for A^'^'^^'^ of (jlSh in section 
12. 3[ The representation ()A.20p , ()A.2ip for the amplitude then also allows the treatment 
of the final state photon 7(5', /u) in the same way as the initial state photon, that is 
replacing it by a dipole in the way explained in section [2j 
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B Decomposition using spin sums 

The free quark propagator for mass niq can be written as follows 



+ '(^°-^°)(i)3 / S^^'^^~^H-^ + -.)l3, (B.l) 



where = (k^ + m^)^/^, and I3 is the unit matrix in colour space. We define the 
Dirac times colour spinors as 

Ur{k) = ux{k)eA, (B.2) 

vr{k) = vx{k)eA. (B.3) 

Here r = (A, A) are the spin times colour indices with A G {1/2, —1/2} and A £ {1, 2, 3}. 
ux{k) and vx{k) are the usual Dirac spinors for mass niq, and are orthonormal basis 
vectors in colour space. We then have for the spin and colour sums 

^Ur{k)Ur{k) = (^ + mg)l3, (B.4) 
r 

Y,Mk)vr{k) = (^-m,)l3. (B.5) 

r 

Inserting this and the Fourier representations for the ^-functions in (jB.ip we get 

\e-'^-''ur{k)ur{k)e'^-y - e^^-^w,(A:)?;,(A;)e-^^-^} . (B.6) 

r 

Here 

k = 



k 



K = k ) ■ ^^-^^ 

We find it convenient to introduce matrix notation. For this we define 

{x\k^,r) = Urik)e~'^'^'' , 

{k^,r\x) = e'^^^'urik), 

{x\k^,r) = Vr{k)e'^'^'' , 

{k^\x) = e-'^^'^Vrik) . (B.8) 

We also use the summation convention for spin times colour indices and momenta, 
implying integration over two identical momenta k with the usual invariant phase space 
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measure. Note that this measure contains the energy denominator and not k^. To 
give an example we write 



d^k 



(27r)32A;0 

In this way we can represent the free propagator (jB.ip as 



Ur{k)ur{k)e'^^y . 



(B.9) 



1 did 
27r uj + ie 



[{x\k^,r){k^,r\y) - {x\ku;,r)(kuj,r\y)] .(B.IO) 



In (jB.lOp the free propagator Sp'^^ is written as the sum of dyadic products of quark 
and antiquark spinors where the energies in the plane wave factors in (jB.SP are off the 
energy shell. 

In the remaining part of this appendix we give the derivation of (j24p . Integration 
by parts and insertion of (jB.lOp for the free propagators leads from ()23p to 



d'^x d'^z d'^z 



Tr 



(B.ll) 



1 



(2vr)2 
Tr 



duj 
uj + ie 



^ du)' 
^uj' + ie 



d^x d^z d^z' 



i^,, +m,)Sf{z',Q;G)^^'Sf{Q,z-G){i^, +m. 



[{z\k^,r){k^,r\x) - {z\k^,r){k^,r\x)\e "'^7'' 

[{x\k'^,y){k'^,y\z') - (x|fc>^)(F;;;7|z')]} , 

where in analogy to ()B.7p we have defined 



(B.12) 



k' 

k.j 



k' ' ' 



We define 

(/|i?('?)^|z) = +m,)sf{z',Q-G)^^sf{Q,z-G){i^, +m,) 



(B.13) 



(B.14) 
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and further 

(fe:,„r'|i?('?)'^|A;^,r) = j z z' {k^, ,r'\z'){z'\B^'^'>^'\z){z\k^,r) , (B.15) 

{ki,,r'\B^'^^>^\k-^) = Jd'zd'z' {ki,,ry){z'\B^'^^'^\z){z\K;;T), (B.16) 

(fc^|i?(^)nA;.,r) = I d^z d^z' {V;;y\z'){z'\B^^^^^\z)^^^ (B.17) 

= J d^z d^z' (^~P\z'){z'\B^''^^'\z)iz\l^ (B.18) 

and 



(A;,,r|i?'^(g)|Cr') = j d'^x {k^,r\x)e-'^'^Y {x\kl, ,r') 

= {2TTf5^*\k^-k'^,-q)Ur{k)YUr:{k'), (B.19) 



{k^,r\E%q)\k'^,,r') = j d'^x {k^,r\x)e-'^^Y {x\k'^, ,r') 

= {2TTf5^*\k^ + kl,-q)Ur{k)YVr>{k'), (B.20) 

{k~^\E%q)\k'^,y) = j d*x{k~^\x)e-'^'-'Y{x\K'y) 

= {27r)^6^^^{-k^ - C - q)vr{kh''ur>{k') , (B.21) 



ik^,r\E''{q)\k'^,,r') = J d^x {k^,r\x)e-'^^Y {x\kl, y) 

= {2TT)^6^^H-k^ + k'^, -q)vr(k)Yvr'{k'). (B.22) 

Note that the matrix elements of E'^ depend on q, in contrast to those of B^'^^^^. We 
find from ([Bl2]) 



A('^^^"'{q) 



duj f duj' 



(27r)2 J Lo + ie J uj' + ie 

- {V;y\B^'^^'' \k^,r){k^,r\E''iq) jF^^) 

+ {kl,,r'\B'^'^^^^\k^,r){k^,r\E''{q)\kl,,r') 
-{k'^,y\B^'^'^^\l—r){K:rr\E''{q)\k'^,y) 

+ (K7?\B^'i)^\k~?){l~^\E% \ . (B.23) 



Inserting the exphcit expressions ([RT9]) to ([R22]) for the matrix elements of E^iq) and 
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restoring all sums and integrations we find 

-6{k° + k'^-q'^+LU + uj')6^^\k + k' - cl){V^\B^''^^'\k^,r)ur{k)Yvr'{k') 
+ -k'^-q^ + LO- uj')5^^\-k - k' - q)(A;^,, r'|5('?)'^|A:^, r)'u,(A:)7''n,,(A;') 

+ <5(-fc° + - g° - w + u;')'^^^^(-k + k' - ci){^^\B^'''^''\k^)vr{k)Yvr'{k')^ ■ 
Performing the cu' and k' integrations leads to 



+ (A:52),r'|S(5)/^|fe^,r)n,(A:)7'^tx,/(A:'(2))(-gO + fc^ - k'^^^^ + uj + ie)'\2k'^^^^)-^ 
- {kf\r'\B^^^^'\k:;;^)vr{kh•'ur>{k'^^''){-q° - k^ - A:'^^)^ -oj + ie)-\2k'^^'>^)-^ 

+ {k'S^^ , r'l ira^r (fc)7"t;.' ) {q^ + A:^ - + + ie)-\2k"^^^'^)-^ 

where due to the delta functions in the uj' integrations the k'^^^ satisfy relations similar 
to (iRTl) and (lBl3|) . Exphcitly, 

k'W =k;(l) = q-k, 



A;'W° = -^m2 + (k'(i))2, (B.26) 

k'(2)=kL(^) = -q+k, 

= ./m2 + (k'(2))2, (B.27) 



kf^' = -q^ + k'+UJ, 
k'(3) =k:(3) = _q_k, 



A;'(3)o = Jm2 + (k'(3))2, (B.28) 



k'W=k:(^) =q + k, 







A;'(4)o = -^m2 + (k'(4))2^ (B.29) 
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C The matrix element Ad^^^ 



In this appendix we discuss the amplitude A4^^^ (jl9p in a similar way as we did for 
M''"'^ in sections \T3\ and [2^ Our aim is to derive (jM]) and ([I5]) . 
We start from (1191) and write Ai^^^ as 



27:mr)Z„ — 

q',q 



y^.Qq'Qq 



p',s'\^\p,s) {-l)Tr j^sl!'\0,0;G) <^"(g)>^, (C.l 



where 



-ipy 



and 



Inserting in (jas]) factors of 1 from 1^ and ^ we get after integrating by parts 



(C.2) 
(C.3) 



d^xd^zd^z'Tv 



Thus A^^^^iq) has the same structure as A^'^')f^^{q) (jB.lip with the replacement 



Sf''\z,x)e-"^^YS^^'''\x,z')\. 



(C.4) 



si?^(z',0;G)7'^4^^(0,z;G) ^ S^pl\z' , z; G) . (C.5) 

The next step is to insert for 5^'^^ in ()C.4p the expansion (|B.10p . Then all remain- 
ing steps done for A^'^'^^^'^ (q) in appendix iBl can be repeated for A^i^^'^{q). With the 
replacement (jC.Sp we define the quantities analogous to (|B.14p - (|B.18p 

(z'|C(«)|z) = i-i^,, +mg)S^^\z',z-G)ii^, +m,), (C.6) 

{ki„r'\G('^')\k^,r) = I d^z d'z'{kl,,ry){z'\C^'^^\z){z\k^,r), (C.7) 
etc. In analogy to (pl]) - (f29]) we get then 



where 



duj 



271 J uj + iej (27r)32/c0 



d^k 



(C.8) 



(C.9) 
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The bj (j = 1, . . . ,4) are defined as the aj in (p6]) - ([29]) with the replacement 



Inserting now (jC.Sp in (jC.ip we get 

4 



where 



2'iTmnZ. 



(C.IO) 



(C.ll) 



(C.12) 
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Explicitly we get with the momenta k'^^^ and k'^^^ as in ()B.26P - ()B.29P 

t' ,r 

Ur{k)Zqj''Vr'{k'^^^) , (C.13) 



G 



q',q 



duj 



(fk 



(2fc'P'")-i J^(7(9',>.),«(t2'>,r'),p(p',<i')|rl«''"|«(t„,r),p(p,s)) 

r',r 



q',q 



d^k 



uj + ie J (27r)32A;0^ ^ ^ 
i2k'(^^^)-'Y.(^iq',f,),qik^,r),qikf\r'),pip\s')\T^^^ 

r',r 

Vr{k)Zg-f''Ur'{k'^^^), (C.15) 



q',q 



doj 



d^k 



(2fc"'"»)-'2(7(9',rt,9(t^.'-),p(p',/)|rl»'''>|{(t2",r'),p(!),s)) 

r',r 

(-l)i;,(fc)Zg7''?;,,(fc'(^)), (C.16) 
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where 



(C.17) 



{l{q\^J),q{k^^\r')Mp\s')\T^'''%{k^,r),p{p,s)) = (C.18) 

\p\s'Mp,s) (-l)TV [7^4^')(0,0;G)1 (fe^^)^ 



G 

fi),qik^,r), q{kf) , r'),p{p', s')\T^'i' \p{p, s)) = (C.19) 
'{p',s'mp,s)TT [7'^4^')(0,0;G)1 (A;53),/|c('')|fc;:;7)\ 



'ZTTTTlp Zp Zq 



G 



M, l^),q{k^,r),pip', s')|T(^' ^) \q{k'S^^ , r'),p{p, s)) = 

{p',s'mp,s)TT\rS^/\0,0;G)\ ik'S^\r'\C^'i^\l~T) 



(C.20) 



I'KTXlpZpZq 



G 



In ()C.13P - ()C.16P we can replace Zq-^^ by (HH). Choosing the momenta in (jC.lSP ac- 
cording to dll]) leads to (|i5|) . 
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